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ABSTRACT 

We study the kinematics of a local sample of stars, located within a cylinder of 500 pc radius centered on the 
Sun, in the RAVE dataset. We find clear asymmetries in the vr-v^ velocity distributions of thin and thick disk 
stars: there are more stars moving radially outwards for low azimuthal velocities and more radially inwards for 
high azimuthal velocities. Such asymmetries have been previously reported for the thin disk as being due to 
the Galactic bar, but this is the first time that the same type of structures are seen in the thick disk. Our findings 
imply that the velocities of thick disk stars should no longer be described by Schwarzschild’s, multivariate 
Gaussian or purely axisymmetric distributions. Furthermore, the nature of previously reported substructures 
in the thick disk needs to be revisited as these could be associated with dynamical resonances rather than to 
accretion events. It is clear that dynamical models of the Galaxy must fit the 3D velocity distributions of the 
disks, rather than the projected ID, if we are to understand the Galaxy fully. 

Subject headings: Galaxy: disk — Galaxy: kinematics and dynamics — Galaxy: structure 


1. INTRODUCTION 

The velocity distribution of stars in the Solar neighborhood 
contains key information about the current dynamical state of 
our Galaxy and also about its history. The kinematics of stars 
can be used to derive both the mass distribution of the Milky 
Way through sophisticated dynamical models, as well as to 
identify accretion events. 

The velocities of thin disk stars are often described by the 
Schwarzschild distribution function, which considers each of 
the velocity components se parat ely. However, data from the 
Hipparcos mission dPerrvman et all 119971 ESA 1997) and 
later from the Geneva-Copenhagen survey (iHolmberg et all 
120091) have revealed a more com plex distribution with signifi¬ 
cant overdensities and stru cture (I Dehnenl 1 998llFamaev et afl 
120051 lAntoia et all 120081) . Some of the over-densities and 
distortions of the velocity distribution appear to be the im- 


prints of the non- 
namely the spira 

axisymmetric components of the Milks 

Way, 

1 arms and Galactic bar (e.g. Dehnen 

2000; 

De Simone et al. 

2004; Sellwood 2010: Antoia et al. 

2011; 

McMillan 2011 

). Streams or moving groups are formed by 


stars on orbits that are close to resonant with the natural fre¬ 
quencies of the spiral arms and/or the bar. Examples of such 
moving groups that are heterogeneous in age and chemical 
composition are the Pleiades, Hyades, Sirius, and Hercules 
streams (the latter very likely driven by the bar). 

* ESA Research Fellow. 


On the other hand, the velocity distribution of the thick 
disk has been studied in less detail thus far because of lim¬ 
itations in the size of volume complete samples. In prac¬ 
tice, this has implied that Gaussian distri butio ns were used to 
fit the kinematics o f thick disk stars (see lBinnev et al.lf2014t 
ISharma et all 120141 for a recent discussion on how Gaus¬ 
sian functions poorly fit all velocity components). Further, 
substructures have also been reported in the thick disk (e.g . 
iGilmore et al.1120021 : iNavarro et all 120041 iHelmi et~aT1120061) , 

identified through statistical comparisons with Galactic mod¬ 
els or with simple kinematic models such as those discussed 
above. 

Many of these substructures have been attributed to accre¬ 
tion events, as these typically leave behind streams of stars 
with similar velocities that do not necessarily appear to be 
spatially coherent near the Sun because of the short mix¬ 
ing timescales in the inner Galaxy. However, recent_model- 
ing has_shown that the impact of spiral arms (ISolwav et all 
1201 2t iFaure et akll2014l) and the Galactic bar on the kinemat- 
ics of stars in the th ick d isk is non-negl igible. For example, 
iMonari et all (120131) and iMonaril (120141) have found in their 
simulations that there is as much resonant trapping in the thick 
disk as in the thin disk. Another clear signature of the impact 
of the bar in their thick disk simulations is a bimodality in the 
velocity distribution for stars located near the Outer Lindblad 
Resonance, similar to that observed in the thin disk. 

Here we explore whether these features are present in local 
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samples of thick disk stars, especially now that such samples 
have increase d in size by larg e fa ctors, (as in e.g. LA MOST 
and SEGUE, ICui et alJ I 20 I 2 I and lYannv et alJ 120091 respec¬ 
tively). For example the Geneva-Copenhagen survey con¬ 
tained ~ 17,000 stars, while ~ 400,000 stars have now 
been measu red by the RAdial Velocity Experiment (RAVE) 
(ISteinmetz et al.ll2006l) . We report on the analysis of the local 
RAVE dataset, and find indeed clear asymmetries/structures 
in the velocity distributions of local thick disk stars, which 
can be attributed to the resonant interaction with the Galac¬ 
tic bar. In Section [2] we present the dataset and the selection 
of the different populations, in Section[3]the analysis, and we 
conclude in Section[4]with a discussion on the implications of 
our findings. 

2. OBSERVATIONS AND DATA SELECTION 

In this study w e use the RAVE Data Release 4 (DR4) 
(iKordonatis et al.ll2013bl) . The stellar atmospheric parameters 
of the DR 4 are computed using tw o di fferent pipelines, pre - 
sented by IKordonatis et all (1201 ll) and ISiebert et al.1 (1201 ll) . 
The stellar parallaxes that we use were obtain ed through the 
Bayesian distance-finding method of lBinnev et al.l (120141) . 

First we select stars with i) signal-to-noise ratio better than 
20, ii) the firs t mo rphological flag indicating that they are nor¬ 
mal stars dMatiievic et al.ll20l2h . and iii) converged algorithm 
of computation of the physical parameters. From these, we 
further select those in a cylinder with radius of 500 pc centred 
on the Sun’s position. This results in a sample of 162153 
stars with 6D phase-space information, of which 76% are 
dwarf stars and 24% are giants. The DR4 proper motions 
were compi led from several catalogs and here we use UCAC4 
(IZacharias et al .1120131). 

Following iReid et akl (120141) we assume that the Sun is 
at X = -8.34kpc and take a circular velocity at the Sun 
of Vo = 240 km s 1 . For the velocity of the Sun with re¬ 
spect to the Local Standard of Rest we adopt (£/©, V 0 , W©) = 
(10,12,7)kms _1 (ISchonrich et al.l2010l) . The resulting value 
of (Vo + V 0 )/Ro is 30.2 km s~' kpc^ 1 , which is compati¬ 
ble with that fro m the reflex motion of th e Sgr A* 30.2 + 
0.2kms _1 kpc 1 (IReid & Brunthaled 120041) . With these val¬ 
ues, we compute the stars’ cylindrical velocities: vr (posi¬ 
tive towards the Galactic center, in consonance with the usual 
U velocity component) and v,t, (towards the direction of rota¬ 
tion). 

From the selected sample we consider 4 different subsets 
of stars based on their height and their metallicity to maxi¬ 
mize or minimize the number of thin or thick disk stars. The 
properties of each subset and relative thin/thick/halo fractions 
are listed in Table |T] Two of the subsets are located on the 
plane but have metallicities corresponding to thin (1) and to 
thick (3) disk components, respectively. The other two are lo¬ 
cated far from the plane and have intermediate (2) and low (4) 
metallicities and could be both associated with the thick disk. 

For each subset it is important to estimate the fraction of 
stars that could belong to a different population than desired. 
We have derived two different estimates of these fractions 
for each of the samples. The first estimate, which we term 
RAVE-fit, is based on an admittedly simplistic three Gaussian 
population model (old thin disk, thick disk and halo) fit to 
the met allicity distribution to RAVE data bv lKordopatis et ahl 
(l2013ai their tables 1 and 2). We use the fits derived for all 
stars with galactocentric radius between 7.5 kpc and 8.5 kpc, 
to estimate the population fraction for samples (1) and (3). 
Since these samples have an additional constraint, namely 


|Z| < 0.5 kpc, the fractions of thick disk and halo stars are 
probably overestimated. For samples (2) and (4) we use the 
RAVE fits derived for stars 1 < |Z| < 2 kpc in the same ra¬ 
dial range. In this case, since in our samples we consider all 
stars with |Z| > 0.5 kpc, it is likely that the fraction of thin 
and thick disk stars is underestimated, while that of the halo 
is overestimated. In fact, if we assume that the halo has no 
net rotation and that all stars with < 0 belong to the halo, 
we can estimate the fraction of halo stars as twice that of stars 
with < 0. We find this to be of only 3%, 0.7% and 4% for 
samples (2), (3) and (4), respectively, i.e. much smaller than 
the fractions obtained through the RAVE-fit. 

The second estimate (fit 2) of the contamination in our sub¬ 
sets is based on a simple model with two populations (thin 
and thick disk) with specified density and metallicity distri¬ 
butions. We use two exponential disks with vertical scale 
heights of h ZJ hm = 0.3 kpc and h ZJ hick = 0.9 kpc and scale 
lengths of fiR^hm = 2.6 kpc and Iir, thick = 3.6 kpc 2 , and a local 
density normalization of 12%, all as measured by lJuric et al.l 
(12008I) 3 . We also assume Gaussian metallicity distributions 
with means <[M/H] rtm > = -0.1 and <[M/H] (tect > = -0.78 
and dispersions <X[m/h], w „ = 0.2 and cr tM /H] rtict = 0.3 (simi- 
l ar to the intermediate-old thin and thick disk populations of 
iRobin et~aT] 120031 respectively). We estimate the fraction of 
each population by integrating between the given ranges of 
metallicities and heights. For samples (2) and (4) we assume 
a maximum height of 1.5 kpc. 

The population fractions estimated with the two methods 
(RAVE-fit and fit 2) indicate that the contamination of thick 
disk stars in sample (1) is very low. On the other hand, sam¬ 
ples (2), (3) and (4) are clearly dominated by the thick disk as 
desired. 

3. STATISTICAL ANALYSIS 

In Fig. [I] we show the velocity distributions of the different 
samples using scatter plots (left) and a kernel density estima¬ 
tor (right, see caption for details). The velocity distribution of 
the thin disk, subset (1), is not homogeneous and de picts o ver- 
densities and asymmetries, as already reported in lAntoia et al.l 
(120121) for RAVE thin disk stars. We see a clear asymmetry: 
stars with v t! , < 220 km s -1 are shifted to the left part of the 
distribution (vr < 0). Interestingly, this asymmetry is visi¬ 
ble in the thick disk subsets for both scatter and density plots 
shown in the remaining panels. 

To study this in more detail, we use the density field shown 
in the right panels of Fig.|T|to compute the difference between 
regions with positive and negative vr. In practise, if <r + is the 
density in a certain pixel (vr, v ( t,) of the grid (of 2 km s -1 size) 
and cr_ is the number for the symmetric pixel (-Vr, v,/,), we 
compute A = cr + - cr_. This quantity is plotted in Fig. [2] 
Red colors indicate A > 0 (more stars with vr > 0), and blue 
colors A < 0. 

The upper left panel in Fig. [2] corresponding to the thin disk 
sample, clearly shows that the region with v# > 240kms~ 1 
has an excess of stars with Vr > 0, while the contrary is true 
for v ( t, < 240 km s' 1 . This separation is not at constant v,j, for 
all vr (i.e horizontal in the v K -v (tl diagram) but rather varies 
with vr. The other three panels, corresponding to the thick 

2 Th e scale lenghts and density normalization of the disks are uncertain, 
see e.g. IBensbv et all201 11) . 

3 IRobin et~aD <20030 gives a normalization of 27% for the intermediate- 
age to old thin to thick disk. This would yield an even lower thin disk con¬ 
tamination in samples (2) to (4). 
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TABLE 1 

Properties of the different cuts. 




N e z 
kpc 

e dist 2 d 

kpc 

e VR 

kms" 

kms 1 

e Vz 

kms -1 

[M/H] 

dex 

RAVE-fit 

thin thick halo 

fit 2 

thin thick 

1 

[M/H] > -0.1 dex 

|Z| < 0.5 kpc 

47883 0.05 

0.05 

5.0 

4.3 

3.3 

0.10 

96 

4 

0.5 

99.8 

0.2 

2 

[M/H] < -0.45 dex |Z| > 0.5 kpc 

5123 0.38 

0.13 

15.7 

15.1 

5.2 

0.10 

0.2 

88 

11 

13 

87 

3 

[M/H] < -0.7 dex 

|Z| < 0.5 kpc 

21624 0.06 

0.06 

6.6 

5.0 

4.0 

0.12 

0.7 

78 

22 

3 

97 

4 

[M/H] < -0.7 dex 

|Z| > 0.5 kpc 

2939 0.40 

0.13 

17.3 

18.0 

5.4 

0.10 

0.003 

68 

32 

0.7 

99.3 


Note. — The first columns show the cuts, the number of stars N, median errors in vertical position Z, in horizontal distance in the plane from the Sun (listed, 
in the velocity components ( vr , v# and vz). and in the metallicity [M/H], The last 5 columns show the thin disk, thick disk and halo fraction of stars computed 
according to different models (named RAVE-fit and fit 2, see text for caveats with respect to the estimated halo fraction). 
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Fig. 1.— Velocities of stars in the different cuts of the R AVE dataset of Table [7] Left column: scatter plots. Right column: density obtained through the 
Epanechnikov adaptive kernel density estimator method flSilvermanll 19861) with an adaptability exponent of 0.1. The density was estimated in a uniform grid of 
2 km s -1 . The black contours indicate the following levels in units of the maximum density: 0.005, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.995. 
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Fig. 2.— Differences in the density field (depicted in the right panels of Fig.fT} for stars with positive and negative vr for the different subsets. The density was 
normalized to the number of stars in each sample and, therefore, the color bars indicate difference in the number of stars in each pixel of the grid (of 2 km s _1 ). 
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Fig. 3.— Left: Difference in the number of counts between positive and negative vr, A = (N+ — N~) (black curves) and the same number scaled to the total 
number of stars in both bins A = (N+ — N-)/(N+ + AL) (red curves, right axis) as function for the different subsets. Right: Mean vr velocity for bins in for 
the different samples. 


disk subsets, depict a similar asymmetry. However, in these 
cases, the asymmetry is not as sharp and clear as for the thin 
disk set. This is probably due to the larger velocity errors and 
to the lower number of stars, especially for samples (2) and 
(4). Nevertheless, the asymmetry is very clear for sample (3), 
and it is located roughly at the same velocities as for sample 
(1). Despite limitations for samples (2) and (4), it is still clear 
that red colors dominate the half upper part of the distribution 
while blue colors dominate in the other half 4 . 

The left panel of Fig. [3] is a one-dimensional version of 
Fig. 0 where we plot for each v (/ , the difference between the 
counts for v R > 0 (N+) and those for v R < 0 (AL). We 
use the m ethod presented in lScarg le et all (120131) and imple¬ 
mented in iVanderPlas et akl (120121) to bin the data in It 
is a non-parametric technique that finds the optimal data seg¬ 
ments of variable size that maximize a certain fitness function 
in a Bayesian likelihood framework and based on Poissonian 
statistics. Although the binning choice is arbitrary, we have 
checked that our conclusions do not change if we use bins of 
equal size or bins with equal number of stars. For this figure 
we plot also in red (right vertical axis) the relative asymmetry 
in the counts, i.e. A = ( N+ - AL)/(A+ + AL). 

4 Note that a wrong assumption of the peculiar velocity of the Sun t/ Q 
would also produce an asymmetry in the counts of vr > 0 with respect to 
vr < 0, but this would be positive or negative everywhere and would not 
depend on as we see here. Note also that a different assumption of the 
values for Rq, V© and Vo can not produce the observed asymmetries, only 
shifting the positions or the velocity . 


This plot shows the trend already highlighted in the two- 
dimensional plots of Fig. [2] There is a large asymmetry in the 
counts towards v« > 0 that peaks at ~ 230 - 250 kms -1 
and extends from v (!> ~ 200 to ~ 275 km s -1 . A smaller but 
significant (note the small errors) asymmetry A < 0 is de¬ 
tected for velocities below ~ 200 km s -1 and at least down to 
120 - 140 kms -1 . This asymmetry is in a region of the ve¬ 
locity plane with a lower number of stars, and is thus better 
seen for the red curves which corresponds to the normalized 
counts A. 

On the right panel of Fig. [3] we show the mean v R as func¬ 
tion of v ( t, for all the subsets using the same binning as before. 
The shape of these curves are similar to those on the left plot. 
This is because an excess of counts for positive v R reflects 
in a positive mean v R , and conversely. The change in sign 
of the mean v R is significant for all four samples and occurs 
at i '0 ~ 200 - 220kms -1 . It is noteworthy that the nega¬ 
tive asymmetry seen both in the counts and in the mean v R 
extends well into low v^, a region more often thought to be 
characteristic of accretion events. For instance, thi s is where 
the Arcturus stream is also located dNavarro et al.ll2004h . 

For the thin disk sample (top panels of Fig. (3} we ob¬ 
serve also other smaller bumps and more detailed features. 
In particular, two clear positive peaks are seen both in the 
counts A and in the mean Vr. These are due to known kine¬ 
matic over-densities in the Solar Neighborhood. Some hints 
of similar bumps (e.g. at ~ 215 - 220 km s ’ 1 and at 
v (!l ~ 250 - 260 km s -1 ) are also seen for the thick disk sample 
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(3). 

The results presented above are robust to the specific choice 
of proper motion catalog, as we have checked that no signif¬ 
icant differences are seen when using other catalogs, namely 
UCAC2 and UC AC3 dZacharias et all 120101 and references 
therei n), Tvcho-2 (iHOuet al J2000h and PPMXL dRoeser et al.1 
120101) . We also obtained the same results when we scaled 
the distances by factors of 0.8 and 1.2, i.e. assuming that the 
distances were overestimated and underestimated by 20%, re¬ 
spectively. Furthermore the choice of the volume of the cylin¬ 
der has no effect on our results: stars with the same cuts but 
located in a cylindrical ring of radius between 0.5 and 1. kpc 
(dominated by giant stars instead of by dwarfs as in our initial 
cylinder) show the same asymmetries. 

As mentioned in SectionQ] the main effect of the bar on the 
local velocity distribution is a bimodality in the velocity dis¬ 
tribution near the Outer Lindblad Resonance (iKalnai si 1199 lL 
lDehnenll2000l) . This produces the Hercules stream, an excess 
of stars with negative Vr at velocities around = 190 km s 1 
(heliocentric velocity V ~ -50 km s _1 ), and the dominant 
mode of low-velocity stars centered around the local standard 
of rest (v^ = 240 km s 1 ) and with an elongation through pos¬ 
itive v R . We believe that the asymmetries that we observe have 
the same origin. 

4. DISCUSSION AND CONCLUSIONS 

We have found clear asymmetries in the vr-v^, velocity dis¬ 
tributions of thin and thick disk stars near the Sun. In the thin 
disk such asymmetries are due to well-known over-densities 
such as the Hercules stream, which has been explained by the 
resonant effects of the bar near its OLR. 

This is the first time that the same type of structures and 
asymmetries are seen in the thick disk. The features are sig¬ 
nificant for the three different thick disk samples considered 
based on metallicity and height above the plane, which we 
have estimated to have low contamination from thin disk and 
halo stars. These findings suggest that the Galactic bar leaves 
strong imprints on both the thin and thick disk. 

It is clear that the observed v ltl velocities are highly skewed, 
not following a Gaussian distribution. iBinnev et al.1 (120141) 
also showed that the radial and vertical velocities can not be 
fitted by Gaussian functions not only due to moving groups 
but also because they peak more sharply than Gaussians. Our 
results also imply that the velocities of thick disk stars should 
no longer be described by independent Schwarzschild, mul¬ 
tivariate Gaussian or purely axisymmetric distributions: the 
reported asymmetry is both v t!l and vr dependent. The RAVE 
d ataset has_also shown peculiar vertical velocity patterns 
(iWilliams et al.1120131) and other studies based on simulations 
have found that the non-axisymmetries of the G alaxy can also 
influence the vertical velocity distribution (iFaure et al.1120 141) . 
It follows that that dynamical models of the Galaxy must fit 
the 3D velocity distributions, rather than the projected ID. 

Simulations show that the deformations and over-densities 
of the velocity distribution caused by the bar change with po¬ 
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sition (both in radius and azimu thal angle) in the thin (e.g. 
lpehnen|2000y Antoia et al.l2014l) and in the thick disk as well 
dMonari et al.ll2013h . This implies that the asymmetries in the 
velocity distribution of the thick disk must change with spa¬ 
tial position. Interestingly, rotatio nal lags and asymmetries 
in t he thick disk we r e repo rted by iHumphrevs et al.l (1201 ll) 
and iJavaraman et al.l (120131) . which may be further evidence 
of this. Future modeling of the velocity distributions must 
therefore be position dependent. One should notice that the 
asymmetry projects differently on line of sight and transverse 
velocity, creating different signatures. 

This also implies that the nature of previously found sub¬ 
structures in the thick disk needs to be revisited as these could 
be associated with dynamical resonances rather than to accre¬ 
tion events. Specifically, the Arcturus stream would seem to 
be naturally explained in this way, being an extension towards 
lower Vc/, velocities, which would also be favored given its 
chemical abundances dWilliams et al.l[2009l) , and hence there 
would be no need to invoke an y acc retion events nor ringing 
due to a past merger event (|Minche y et al.150 09). Also the 
substructure reported by iGilmore et al.l (l2002i) could perhaps 
be explained along similar lines, although it was su ggested 
that th is is part of a metal-weak thick disk (Kordonat is et al.l 
12013d) . The role of the bar on the formation of such structures 
should thus be investigated. 

It is clear that with the advent of larger and more precise 
samples of the disk kinematics we are entering a new era 
where the classic velocity distribution models are not suffi¬ 
cient and the effects of the non-axisymmetry in the disk have 
to be taken into account in the modeling. This is particularly 
relevant now given that in approximately two years time the 
first Gaia data will be published and we expect to detect these 
and perhaps other asymmetries far beyond the Sun. 
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